Field Theoretic Description of High Energy 
Neutrino Interactions in Matter 



G. Domokos and S. Kovesi-Domokos 
The Henry A. Rowland Department of Physics and Astronomy 
The Johns Hopkins University 
Baltimore, MD 21218Q 

a^ : 
a^ . 

', In this paper we begin the development of a formaUsm for the description of high energy 

Q I neutrino interactions. It is based upon field theory quantized on a null plane. We set up 

Q ' the general formalism as well as some techniques needed to perform phenomenological 

■ calculations. We show that the formalism developed by Wolfenstein is recovered at 

CO . the cost of making two approximations: one has to treat the charged lepton fields in 

I , the Hartree-Fock approximation and one has to take the short distance limit of the 

I Hartree-Fock correlation function. As an example, we discuss the resonant interaction 

' of electron neutrinos in an electron gas. 

cn ■ 

(N : 

11^ 1 Introduction 

'NT ■ 

a^ ; 

^ I The theory of neutrino interactions at high energies will play an increasingly 

Q^l important role in the future. High energy neutrino oscillation and interaction 

p • experiments will be performed by means of accelerator generated neutrino beams 

^ ■ as well as in experiments relying upon extraterrestrial sources of neutrinos such 

• • . as active galactic nuclei (AGN) and binary systems of stars. In the latter setup, 

I the center of mass energies in a neutrino interaction with matter are likely to 

^ ' reach a few TeV, corresponding to laboratory energies of the incident neutrinos 

' of the order of a few PeV. Traditionally, neutrino interactions have been either 

described by the single gauge boson exchange approximation (as in the case of 
the theory of deep inelastic scattering of neutrinos on quarks) or in a contact 
interaction approximation (for instance, in the theory of neutrino oscillations in 
matter or in vacuum as developed by Wolfenstein, [jl|.) While there exist several 
derivations of the Wolfenstein formalism 0, |^, it is clear that a more general 
formalism is needed at the highest energies. To quote but one example, if electron 
antineutrinos of an energy of about 6.4 PeV pass through an electron gas, they 
excite the W~resonance (the "Glashow resonance"). Such a situation may oc- 
cur, for instance, when electron antineutrinos generated deep in an AGN pass the 
electron plasma on their way out. Clearly, neither the contact interaction approx- 
imation, nor a gauge boson exchange approximation are adequate; the resonance 
is excited in the s-channel; in this case, one has to be able to treat the presence 
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of matter adequately. Due to the fact that neutrinos are nearly massless and we 
are interested in their interactions at high energies, the null plane or front form 
formulation of the theory is the most convenient one for the purpose. In the next 
section we briefly review the formalism. Thereafter, we outline the procedure for 
obtaining an effective theory of neutrino interactions in matter; as an application, 
we show how Wolfenstein's equation is recovered form the general theory. In sec- 
tion 4 we illustrate the use of the theory by describing the interaction of electron 
antineutrinos (V^) with an electron gas in the energy region where the W boson 
is excited as an s channel resonance. Sec. 5 contains a discussion of the results. 

2 A Review of the Null Plane Formalism 

The null plane or front form of quantum field theory is, in essence, a constrained 
Hamiltonian formulation of a field theory given by its Lagrangian and either oper- 
ator quantization rules or path integral prescription. Instead of a spacelike surface, 
however, as it is the practice in setting up a traditional Hamiltonian formalism, 
one prescribes Cauchy data on a plane with a null normal vector. Such planes al- 
ways contain characteristic lines of a relativistic wave equation. As a consequence, 
the number of independent Cauchy data is smaller than the ones one can prescribe 
on a spacelike surface, see, e.g. This is expressed in the form of constraints 
obeyed by the fields entering the theory. As in any Hamiltonian formulation of a 
relativistically invariant field theory, manifest Lorentz invariance is lost. In what 
follows, we describe the formulation of the theory of a Dirac fermion interacting 
with an external gauge field. The Dirac field may carry a representation of an 
internal symmetry group. The generalization to the case of a Yang-Mills theory 
or a theory containing scalar fields as well is straightforward and it has been de- 
scribed in a number of articles on the subject; a sampling of some articles is given 
in ref . ||^ . Recently, an attempt has been made to restore some of the symmetries 
(parity invariance in particular) in a null plane formulation, see 0. Due to the 
fact that we are interested in a parity violating theory, those developments are 
not needed here. 

We begin with introducing a coordinate transformation in Minkowski space. 
We always work with a metric chosen as g^,^ = diag{l, —1, —1, —1) in the usual 
Cartesian basis. Now introduce the coordinates. 




Correspondingly, the components of the metric tensor become, 

gzt = gtz = 1, gAB = -Sab (2) 

and all other components vanish. Here and in the following, capital Latin sub- 
scripts and superscripts refer to the directions perpendicular to the chosen con- 
jugate null directions given by eq. (|l|). Correspondingly, the generators of the 
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Dirac algebra are, 7* = (7° + 7^) /a/2, 7^ = (7'^ — 7'^) /a/2 and 7"^. The Dirac 
Lagrangian is of the usual form, 

>C = '0 Q<9 - - ^/^ + ?/'m?/' (3) 

Here, d := 1/27'^ 9/i,^± := 7'^ (1 ± 75) ^j^^-*- Internal symmetry indices are 
suppressed. 

We now introduce spinor projectors corresponding to the two null directions 
given in eq. (|1|), viz. 

Pt = 7t7* and P, = 7^7^ (4) 
For the sake of brevity, we also introduce the notation, 

Ptip = and Pzip = X- 
A straightforward manipulation leads to the useful relation: 

75 = -nS' {Pt - Pz) ■ (5) 

Consequently, the chiral projectors, (1 ± 75) /2 act as helicity projectors upon the 
spinors projected to the conjugate null directions given by eq. (|I]). We define the 
helicity projectors, 

^ iif^ (6) 

With this, eq. @) becomes: 

a/2£ = i(j)^dt(j) + ix'dzX + 4>^Y^X + X^1^^4> 

+ xV7^Q5a-^^?(+)AJ^+-^^7(-)Aa^-)0 (7) 

It is to be emphasized that eq. (|^) is completely symmetric under the inter- 
change of the two conjugate null directions. However, the symmetry is destroyed 
if we decide to solve an initial value problem by specifying the Cauchy data on 
one of the null planes. For the sake of definiteness, we describe the procedure by 
specifying initial data on a plane t = 0; the procedure for Cauchyi data specified 
on 2; = is completely analogous to the one described here and it can be obtained 
by interchanging (f) with x- 

We recognize that if t is regarded as the "time" variable, only cf) obeys an 
equation of motion; there is no time derivative in the equation obeyed by x- 
Consequently, the equation obeyed by x is a constraint and thus x can be elimi- 
nated altogether from the equations of motion. Given the the Lagrangian, eq. (|^), 
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the constraint can be solved at least formally. The resulting equations of motion 
obeyed by (f) are, in general, non-local; nevertheless, they are legitimate dynamical 
equations. 

In the case of eq. (^, we proceed by choosing the gauge, = 0. In this gauge 
the solution of the constraint is trivial, since one just has to invert the operator 
dz- In order to present the result, we introduce the covariant derivatives, 

ik = A,t) (8) 
After elimination of the constraints, the end result reads: 
I/V2C = (t)Hvt^H-(f) + (j)H\/i'^^H+(l) 

- ((m - zV(+)^) H+4>) iV^^'^^) H+4> (9) 

In this equation, we denoted k = idz and, indeed, in the gauge chosen the 
easiest way of eliminating the constraint is by means of a Fourier transformation 
in z. The singularity at = has to be eliminated by taking the pricipal value, 
due to Hermiticity requirements. 

We notice that eq. (|) is of the canonical form, 

C = ndt(j) - n, (10) 

as described recently by Jackiw p. In fact, the entire procedure of elimination of 
the constraints follows the pattern described in that reference. 

Were there no mass terms and no mixing between left and right-handed com- 
ponents, there would be a one to one correspondence between particles (antipar- 
ticles) and negative (positive) helicities, respectively. Therefore a two component 
theory as described by eq. (||) would be exact. However, at sufficiently high en- 
ergies one expects that the amplitudes of the "wrong" helicity components are 
suppressed by factors of the 0{ma/E), where ma is some of the eigenvalues of the 
mass matrix appearing in the previous equations, for instance, in eq. (^). There- 
fore, the two component theory is expected to be a good effective theory. Parity 
conservation (see ref. [0) is not an issue since we are dealing with electroweak 
interactions. (We recall that under parity, t ^ and (p ^ X- Thus, 

specifying Caucy data on a t = const, {z = const., resp.) surface, parity is, by 
necessity, violated.) 

The Feynman rules can be easily read off from eq. (^. One notices that 
they have the appearence of the rules for a nonrelativistic Schrodinger theory in 
2+1 dimensions. This is not an accident: it is due to the fact that the stability 
group of a null direction in Minkowski space is E{2), the Euclidean group in two 
dimensions. 
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3 Neutrino Interactions in Matter 



In this section we begin the development of the formahsm needed to describe 
neutrino interactions in matter. For the sake of simphcity, we discuss exphcitly 
the case of neutrino interactions in a uniform electron gas. As we proceed, it 
will become obvious that the formalism can be generalized in a straightforward 
manner to describe other physical situations of interest. 

Our starting point is the effective action, see, e.g. ref. 0, F, vieved as a func- 
tional of classical neutrino and electron fields, denoted by ip and ^ , respectively. 
Both Ip and \l/ are regarded as Dirac fields. Internal symmetry indices are sup- 
pressed as before. 

The effective action can be expanded in a functional power series (Volterra 
series). We argue that in many cases of physical interest, the Volterra series of F 
can be broken off after the first few terms. Clearly, the the second derivatives with 
respect to ip and ^ give the free actions for these fields, with the appropriate mass 
terms. Fourth and higher derivatives give the effective interaction kernels. Both 
from a dimensional argument and from the explicit calculation that follows, one 
realizes that terms proportional to (\&^)*^ (with any arrangement of the space- 
time arguments of the fields) are proportional to n^, Uf. being the density of 
the electron gas. As a consequence, in most cases, terms proportional to higher 
powers of the density can be dropped^. On a similar basis we omit derivatives 
higher than second in because we assume that the neutrino beam considered 
is sufficiently dilute so that self interactions of neutrinos can be neglected. (In 
addition, in the electroweak theory terms containing higher powers of the electron 
density also contain higher powers of the fine structure constant and, perhaps, of 
inverse gauge boson masses.) 

With this in mind, we now write the effective action: 

F[^,vl/] = S'o[^] + ^o[^] + /^(Tm2)H(l,2;3,4)v|/(3)^(4) 

+ j ^(l)vI/(2)/C(l,2;3,4)vI/(3)V^(4) + ... (11) 

In eq. (|ll]), the space-time points have been denoted simply by Arabic numer- 
als; the integrations extend over points occurrring twice under the integral sign. 
Clearly, eq. ( pT]) defines an interacting system of electrons and neutrinos. As- 
suming that Ti and K, have been computed in some approximation, one can solve 
this classical field theory in order to represent the neutrino electron interaction. 
The relevant Green functions can be represented as a sum of tree diagrams, with 
interaction vertices given by the kernels Ti and /C. An easy way of generating 
these diagrams is to introduce an auxiliary second quantization of the electron 
and neutrino fields and then compute the tree diagrams in this theory. 

^ This is, in essence, the classical statistical argument for the dominance of binary collisions 
in a gas of moderate density. 
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As an application, let us show how the Wolfenstein formalism is recovered in 
this framework. 

The crucial point is to make a Hartree-Fock approximation in eq. (0). This 
means that one has to replace bilinear products of the form \E'(1)^(2) in the terms 
proportional to Ti and K, (and in the higher order terms) in the effective action 
by their expectation values in the electron gas. In this way, the effective action 
reduces to a quadratic functional. On performing the above-mentioned trunca- 
tion on the term proportional to ?i, one can absorb the resulting functional into 
Sq\^] to a good approximation. In fact, it follows from straightforward invariance 
arguments that the expectation value ( denote it by H) must be of the form, 

H = y ^/lo^ + + d^^h2d^m + . . . , (12) 

where all fields are to be taken at the same point. 

If the electron gas is non-relativistic, as we assume, the higher derivative terms 
may be omitted and we see that the effect of the self interaction within the electron 
gas is just a mass shift and a wave function renormalization. (In standard many 
body theory, see e.g. ref. 0, one often introduces an energy dependent effective 
mass in order to take some of the higher derivative terms into account. While 
this can be done within the framework of the present formalism, it is unimportant 
from the point of view of the argument that follows.) 

Let us now concentrate on the electron-neutrino interaction term. On replac- 
ing the bilinear electron operator by its expectation value, the resulting expression 
is of the form: 



^ V^(1)K(1,2)V^(2), (13) 
where the kernel is given by 

K(l,2) = j /C(1,3;2,4)C(3,4). (14) 

The quantity denoted by C(3,4) is the two point correlation function of the elec- 
tron gas. Due to the homogeneity of the electron gas, it depends only on the 
difference of its two arguments. On making a non-relativistic approximation to 
the kinematics for quantities referring to the electron gas, one obtains the well 
known expression: 

n( +\ ^expiirmjpp) , . /■ 2 / n f,^^ 

C(x,r)=ne / auusmup expnu Tpp/me) (15) 

2p J -I 

In eq. (0), is the Fermi momentum; r and p are the time and radial distance 
measured in units of the Fermi wavelength, i.e. r = tpp, p = rpp. (In the usual 
nonrelativistic treatment, the exponential factor before the integral is absent; this 
is merely a question of the definition of the chemical potential.) Finally, me stands 
for the value of the effective electron mass in the gas. 
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One can now proceed to decompose the Dirac spinor describing the neutrino 
according to its projections onto the the two conjugate null directions as described 
in the previous Section and eliminate the constraint. In order to present the result, 
we use an operator notation, such that, for instance: 



K0:= yK(l, 2)0(2) 
In this way we obtain the effective Lagrangian: 

- <p^-iV -m + K)(id, + -i^Ky^ {i^ -m + K)(t), (16) 

where V = 7^9^. 

In essence, the variation of eq. (0) yields the Wolfenstein equation describing 
the propagation of a neutrino in a medium. In order to obtain the form usually 
quoted in the literature, e.g. in ref. [|lOl, the following steps are needed. 



1. One assumes that the propagation is one dimensional; hence, can be 
chosen to be zero by an appropriate choice of the coordinate system. 

2. The interaction in the medium is given by the standard model of electroweak 
interactions, i.e. in the rest frame of the gas, 

K = ir,7° + ir„7V (17) 



3. One neglects the "wrong helicity" components, i.e. one approximates, 

H~(l) ^ (f), H^cj) ^ 

and {mutatis mutandis) similarly for the conjugate problem in which z is 
regarded the time and x the dynamical variable. 

4. One realizes that in eq. (p!5| ) the characteristic momentum scale is given by 
Pf, while in K it is the gauge boson mass, M. In any environment of interest 
(with the possible exception of the very early Universe) M ^ pp. 

Hence, one can safely approximate, 

C(x,t) ^ C(0,0) =ne 



As a result, considerable simplifications occur. In particular, after neglecting 
the "wrong" helicities and replacing the the electron correlation function by the 
electron density, one realizes that the coefficients i^^ and Ka occur only in the 
combination: 

K = n, {K, - Ka) (18) 
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By eliminating the constraint, one identifies the the effective Hamiltonian for the 
field 0: 



Heff = 0^ 



If the denominator in eq. ([TUp were just equal to k, this would be identical to the 
Hamiltonian for the Wolfenstein equation: assuming one dimensional propagation 
and choosing the coordinate system appropriately, one would obtain the equation 



discussed e.g. in ref. [0. Under certain circumstances, one is indeed justifed in 
approximating k — k k. However, the correct formula is given by eq. (0) and 
care is needed in cases when the denominatior may become small. 

The exercise just described is useful because it makes the limitations of the 
Wolfenstein formalism explicit. Within the framework where it is usually applied, 
i.e. the treatment of the solar neutrino problem, the limitations are totally in- 
significant. However, often one finds applications of the formalism to problems 
involving high energy neutrino interactions, where it is, at best, of limited validity. 
In particular, one notices that the Hartree-Fock approximation one has to make 
in order to reproduce Wolfenstein's results describes the interaction of neutrinos 
with an electron gas entirely in terms of the creation and subsequent annihilation 
of a hole in the Fermi sea. At energies of interest in high energy neutrino interac- 
tions, with ^/s perhaps a few hundred GeV, large momentum transfer processes 
are important; in particular, electrons can be ejected into the continuum instead 
of annihilating with a hole. 



4 Neutrino Interactions at the W Resonance 

As an illustration of the formalism developed in the preceding Section, we con- 
sider z7^ of a laboratory energy approximately equal to 6.4 PeV, incident upon an 
electron gas. In a previous paper, we discussed the physical circumstances 



under which this process leads to somewhat unanticipated results. However, there 
we used the Wolfenstein formalism uncritically. Here we show that the present 
formalism reproduces the the results of ref. . 



We notice that at the energy mentioned, W is excited as an s-channel reso- 
nance. As a consequence, it is reasonable to retain only those diagrams which 
resonate at the W mass. For the sake of simplicity, we also make the customary 
approximations, viz. we evaluate the self energy part of the W propagator and 
the vertices at the the mass of the W. (For most purposes, this is an adequate 
approximation and, in essence, it is equivalent to using a Breit-Wigner formula 
to describe the resonance.) 

In this approximation the effective electron neutrino interaction can be written 
in the form: 

r,„, = / | ^^^l^^(l)A (1 - 2) V<2)7^^^^(2) + {conj.), (20) 
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where is the couphng evaluated at the resonance and A stands for the propa- 
gator of the W: it is the Fourier transform of the quantity, 



M2 - fc2 - zM7 ■ 

As usual, we denoted the self energy part evaluated at resonance by — iM--^. 
Hence, M is the physical mass and 7 is the total width. By performing a Fierz 
transformation, the last equation can be brought to the more convenient form: 



int 



/ v[/(i)^^_^vl/(2) A (1 - 2) V;(2)7'^^-^^(l) + {conj.) (21) 



In this form one can easily make a Hartree-Fock approximation, by replacing the 
bilinear quantity in the electron operators by its expectation value. On assuming 
that the electron correlation function is evaluated in the rest frame of the gas and 
replacing the correlation function by its value at the origin, the contribution is 



of the form as in eq. (|18D , with = —Ka- At this pont, we are now ready to 
write down the effective Hamiltonian with the resonant interaction. Due to the 
fact that the effective action was reduced to a quadratic functional and we con- 
sider one dimensional propagation only, it is convenient to Fourier transform the 
effective action. In this manner, the computation of the inverses of the operators 
entering the constraints becomes trivial. We write the Hamiltonian in the form 
(c/. eq. (0)): 

2{p-K) 

Here, p stands for the variable conjugate to z; we suppressed the argument of 0; 
integration over p is understood. 

The question arises whether this Hamiltonian can be replaced by a conven- 
tional one, i.e. 

H = 0T 0. 

2p 

We now proceed to show that in the case of the resonant interaction discussed 
here, such a replacement is well justified. First of all, it is convenient to rewrite the 
expression of k in terms of the elastic and total widths of the resonance. This is 
easily accomplished by noticing that the usual invariant variable s can be written 
as s = 2pme, Me being the mass of the electron. Straightforward manipulations 
then lead to the expression near s = M^: 

In eq. (^2]), m-c is a mass scale characterizing the electron gas; it is defined by 
the relation = ne/mg. The following Table, taken from ref. |]T2[, contains the 



values of the characteristic mass for some environments of interest. 
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Electron densities and characteristic masses 
for some environments 



Environment 


?7,c[cm ^] 




stellar interior (sun) 
Earth 

water 


1027 
1.6 X 10^4 
3 X 10^3 


2 X 10^ 

3 X 10^ 
5 X 10^ 



We now notice that there are no vanishing denominators in the Hamiltonian, 
since k, is complex. Furthermore, the correction to the conventional Hamiltonian 
arising from eq. (^) is of the order of magnitude 2K,/p. Near resonance this is 



2k 



P 



7em^me2 
47M2 ■ 



Even for an environment like a stellar interior, the corrections are minuscule. 
There may be substantial corrections to the naive form of the Wolfenstein Hamil- 
tonian in very dense environments, e.g. in the interior of a neutron star or in 
the early Universe. However, in all probability, the Hartree-Fock approximation 
breaks down before these corrections become truly significant. 

One can include flavor degrees of freedom without any difficulty and discuss 
neutrino mixing near the resonance as it has been done in ref [12|. None of the 
conclusions about the validity of the Wolfenstein formalism is altered by such a 
generalization and therefore we shall not dwell on this topic any further. 



5 Discussion 

The basic purpose of this work has been to establish a formalism leading to an 
effectice theory of neutrino interactions in matter. It became clear that the use 
of the null plane formalism is suitable for the description of these interactions. 
Clearly, however, what is needed is a reliable calculation of the effective action. 
There are several possible approaches to this problem with an increasing degree 
of complexity. In the approach discussed here, one puts the classical ffelds corre- 
sponding to all but a few species equal to zero. (Here we discussed the simplest 
case: only the fields corresponding to neutrinos and the target medium are kept.) 
Even at this level, some interesting generalizations are possible; in particular, at 
high neutrino densities, such as in the early Universe, one can no longer neglect the 
interaction between the various fiavors of neutrinos; consequently, one has to deal 
with a rather non-trivial multichannel problem. In some recent articles, Kost- 
elecky, Pantaleone and Samuel, ref. |11| addressed this problem entirely within 



the framework of the Wolfenstein formalism. However, those authors, in essence, 
make a dilute gas approximation by neglecting the effects of the Fermi sea, while 
keeping the interaction terms between neutrinos. The internal consistency of such 
a procedure is not quite obvious and further studies are needed in order to clarify 
the issues involved. 
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More importantly, in discussing high energy neutrino interactions, one has to 
go beyond the Hartree-Fock approximation as already discussed, since in that 
approximation a substantial amount of physical information is lost. Furthermore, 
in order to obtain a rehable theory, one will have to follow the evolution of several 
components of the system even if ultimately only the transition probabilities for 
the stable particles are kept. 
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